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1. Introduction 


An understanding of the structure and evolution of atmospheric 
circulation is ultimately linked with knowledge of atmospheric heat sources 
and sinks. Prior to the Global Weather Experiment (GWE) , from December 1978 
through November 1979, detailed global analyses of heating were lacking due to 
the paucity of data over much of the earth. The enhanced observational system 
of the GWE in combination with global data assimilation models developed for 
the GWE has provided unprecedented global data sets for the study of the 
planetary circulation (National Research Council, 1985). A determination of 
the global distribution of heating during the annual cycle of the GWE is now 
feasible . 

The large-scale distributions of heating have been estimated by several 
methods. In one method, net heating is determined from independent estimates 
of the individual components associated with radiation, condensation and. 
sensible heating (e.g., Budyko, 1963; Newell et al . , 1969, 1974; Otto-Bliesner 
and Johnson, 1982; Zillman and Johnson, 1985). Net heating has also been 
obtained from analysis of the energy balance required by the First Law of 
Thermodynamics (e.g. , Hantel and Baader , 1978, Geller and Avery, 1978, Lau, 
1979) and from the quasi -geostrophic balance of vorticity and energy (e.g., 
Wiin-Nielsen and Brown, 1962; Brown, 1964). 

Numerous investigations have estimated global heating distributions with 
data from the GWE. Kasahara and Mizzi (1985), Masuda (1984), Sardeshmukh 
(1984), Holopainen and Fortelius (1986), Chen and Baker (1986) and Kasahara, 
Mizzi and Mohanty (1987) employed the First Law of Thermodynamics. Boer 
(1986) and Johnson (1985) estimated vertically integrated heating from total 
energy balance. Wei, Johnson and Townsend (1983), Johnson and Wei (1985), 
Johnson, Yanai and Schaack (1987) and Schaack, Johnson and Wei (1990) computed 
heating through vertical integration of the isentropic equation of mass 
continuity. Although exact comparisons of heating distributions from the 
above investigations have not been carried out, overall agreement is apparent 
in the general features of the distributions calculated from GWE data sets. 

The purpose of this atlas is to provide estimates of the global 
distribution of atmospheric heating throughout the annual cycle of the Global 
Weather Experiment. Heating rates are calculated through integration of the 
isentropic mass continuity equation (Wei et al . , 1983; Schaack et al . , ) 

utilizing the European Centre for Medium-range Weather Forecasts ( ECMWF ) Level 
Illb data set. Distributions of monthly, seasonally and annually averaged 
heating are presented for isentropic and isobaric layers within the 
troposphere and for the troposphere as a whole. The horizontal distributions 
of heating identify planetary heat sources and sinks and their temporal 
evolution. A systematic vertical variation of heating that is linked with the 
distribution of continents and oceans is also evident (Schaack et al., 1990). 
Meridional cross sections of temporally and zonally averaged heating are also 
presented . 
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2 . Da ta 


The data used in the preparation of this atlas are from the ECMWF GWE 
Main Level Illb data set for the period December 1978 through November 1979. 
Here "Main" denotes the original analysis of the GWE data. ECMWF has 
subsequently produced a reanalyzed data set for a portion of the GWE. The 
Level Illb data, resulting from assimilation, represent a combination of 
observed and model generated information. The characteristics of the ECMWF 
data assimilation system utilized in processing the GWE data are discussed by 
Bengtsson et al . (1982) and Bengtsson (1983). The characteristics of the 
various assimilation systems employed in the preparation of GWE Level III data 
sets are summarized in tabular form by Daley et al . (1985) and Daley (1987). 

Global distributions of zonal and meridional wind components, temperature 
and geopotential height were available for 0000, 0600, 1200 and 1800 GMT of 
each day of the GWE year on a 1.875° latitude -longitude grid for 15 mandatory 
pressure levels between 1000 and 10 mb. From these data, global isentropic 
analyses were generated on a 3.75° latitude -longitude grid for 0000 and 1200 
GMT of each day of the GWE year through linear interpolation of the relevant 
meteorological variables with the R over Cp power of pressure (R is the gas 
constant for dry air and c„ the specific heat of dry air at constant 
pressure). In employing the Lorenz convention (Dutton and Johnson, 1967) on 
underground isentropes ( i . e . , 9 < 0 S ; 6 S is the surface potential 
temperature), the pressure and wind components were set equal to their 
respective surface values. 

Information at the earth's surface was not provided in the ECMWF Main 
Level Illb analyses. Therefore, surface values of pressure, temperature and 
zonal and meridional wind components were estimated from GWE Level III 
isobaric data. Surface pressure was interpolated by assuming that 
geopotential varies quadratically with the logarithm of pressure, ln(p) 
(Gerrity, 1977). Surface temperature and horizontal wind components were 
assumed to vary linearly with ln(p) except when surface pressure was greater 
than 1000 mb, in which case temperature was extrapolated downward from 1000 mb 
using the standard atmosphere lapse rate of 6.5 K km' 1 . 

The fields of zonal and meridional wind components and geopotential 
height in the ECMWF GWE Main Level Illb data set were archived before 
initialization, while the temperature fields were archived after 
initialization (Bengtsson, 1983). Because of a preference to use only 
uninitialized information in the analyses, temperature was computed from 
thickness and layer mean relative humidity. Due to the possibility of errors 
in the height fields at 150 mb (Julian, 1983), the height information at this 
level was not utilized in creation of the isentropic data set. 

A varying vertical resolution in the isentropic domain was adopted; 10 K 
between 220 K and 380 K, 20 K between 380 K and 500 K followed by 550 K and 
650 K. After a preliminary study of the heating distributions for the months 
of January and July 1979, the 400 K isentropic level was designated as the top 
level (0 T ) for the analysis of heating. The global distributions of heating 
for July calculated using 400 K and 650 K as the top of the isentropic data 
are very similar. However, the inclusion of the data above 400 K appeared to 
have a negative impact on the global distribution for January. 
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Distributions of heating in layers above 400 K in January (not shown) 
which for the most part are located in the stratosphere, appear o con 
unrealistic features in some regions, particularly over eastern A ® la , 
results from large values of divergence in stratospheric isentropic layers in 
some regions which are thought to be unrepresentative (National Research 
Council 1985; Johnson, 1986). Several changes affecting stratospheric 
analyses were made in the assimilation system at ECMWF during the analysis of 
December and January of the GWE year which likely impacted the results for 
isentropic layers above 400 K. The time-averaged pressure of the 400 K 
isentropic level is less than 160 mb over the globe in each month of the GWE 
yea? In January the time-averaged pressure of the 400 K isentrope slopes 
from near 85 mb in the tropics to 150 mb in extratropical latitudes. 


3. Description of Analyses 

The analyses of the annually, seasonally and monthly averaged heating and 
pressure for the GWE year presented in this atlas consist of t e o owing. 

1. Global analyses of vertically averaged heating for the surface to 
400 K isentropic layer. 


2. Global analyses of isobar ic layer -averaged heating, 

(A) surface to 800 mb (C) 600 to 400 mb 

(B) 800 to 600 mb 0 » 400 to 200 mb 

3. Global analyses of isentropic layer -averaged heating; 

(A) surface to 290 K (C) 310 to 330 K 

(B) 290 to 310 K <°) 330 to 400 K 

4. Meridional cross sections of isentropically and isobarically zonally 
averaged heating; meridional profile of zonally and vertical y 
averaged heating. 

5. Global analyses of surface pressure and the pressure difference 
between the earth's surface and the 400-K isentropic surface. 

6. Global analyses of the pressure difference between the upper and 
lower surfaces of the isentropic layers in (3). 


The following points concerning the analyses are noted: 

l xhe heating distributions were calculated from time-averaged data 
through vertical integration of the isentropic mass continuity 
equation. See Appendix A for a detailed discussion of the metl hod o 
calculation and the following subsection for the definition of the 
mass -weighted averages. 
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2. The^relation between the heating rate as presented in this atlas 
b y ® snd the rate of specific heat addition, Q, is given by 

e - q / n 


where n, the Exner function, is equal to c p (p/p 00 ) K with K-R/c p . 

3. All horizontal heating distributions have been filtered to emphasize 
wavelengths greater than 4000 km. Details of the filtering employed 
are discussed in a subsequent subsection along with smoothing of the 
zonally averaged heating distributions. 

4. In the analyses, contours are omitted where the atmospheric layer 
does not exist (i.e., there is zero mass in the layer). The quasi- 
horizontal line at the bottom of the analyses in the zonally averaged 
isentropic cross sections represents the coldest surface potential 
temperature occurring along the latitude circle during the time 
period considered. 

Since the heating was calculated by integration of the isentropic mass 
continuity equation, the most direct method to display the vertical 
distribution from these results is in isentropic coordinates. However, for 
utilization of the results in studies employing isobaric coordinates, 
vertically averaged estimates are also presented for four isobaric layers. 
These distributions were obtained by interpolating the time -averaged 
isentropic profiles of heating to isobaric levels at 50-mb increments between 
200 mb and the earth's surface assuming a linear variation with respect to the 
R over c p power of the temporally averaged pressure. These interpolated 
values were subsequently mass -weighted vertically averaged over the four 
isobaric layers: surface-800 mb, 800-600 mb. 600-400 mb and 400-200 mb. The 
global distributions are expressed as heating per unit mass in order to 
properly account for the condition that the earth's surface bounds the lowest 
layer, and also occasionally bounds the second lowest layer in regions of 
elevated orography. 

Along with the distributions of heating, several additional fields are 
provided. These are the time-averaged surface pressure, the time-averaged 
pressure difference between the surface and the 400-K isentropic surface and 
the time -averaged pressure difference between the upper and lower surfaces of 
the four isentropic layers. These analyses indicate the time-averaged mass 
distribution within the various layers for which heating is estimated. The 
time-averaged pressure of an isentropic surface bounding any of the layers may 
be determined from the analyses by subtracting the sum of the pressure 
*^^^^ e ^ ences the individual isentropic layers up to the level of interest 
from the time-averaged surface pressure. 
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3.1. Definition of Mass -weighted Averages 

With the hydrostatic assumption, the mass within an isentropic layer is 
determined by the pressure difference between the upper and lower isentropic 
levels through integration of the hydrostatic relation 

pj e - -g' 1 dp/dQ, ^ 

where p is density, J 0 -|3z/3e|, p is pressure, 6 is potential temperature, z 
is geopotential height and g is the acceleration due to gravity. The mass- 
weighted temporally and vertically averaged heating is defined by 


"e.t 


0 - ; eu j t 2 pJe0 dtde / r u rve 

0T tl ®L 


e U r t2 


dtd0 


(3) 


The upper and lower isentropic surfaces bounding the isentropic layer are ©u 
and e L , while the ending and beginning of the time averaging period are t, and 
ti, respectively. In the vertical average over the entire extent of the 
isentropic analyses, 0y — ” ^00 ^ and ®L ” ®s ^ ^ ' 

The mass -weighted temporally and zonally averaged heating is defined by 


/A , 

* 2 ti t? • 2n t 2 / / \ 

0 - f f p J a 0 dtdA / J J pJe dtdA , 

0 ti 0 ti 

while the temporally, zonally and vertically averaged heating is defined by 
A 8 At 

e _ J® T J^J%j e 0 dtdxd© / J^ T /Q n /^P J e dtdxde • (5) 

The contribution to the above integrals from underground isentropic layers 
(i.e., where 0 < 0 S ) is identically zero due to the hydrostatic assumption an 
the use of the Lorenz convention (Dutton and Johnson, 1967). 

3.2. Filtering of Heating Distributions 

All horizontal distributions of heating presented in this atlas have been 
filtered in order to isolate large-scale features and suppress noise. First, 
zonal harmonics with wavelengths less than 4000 km are removed through 
truncation of one dimensional Fourier transforms (Wei et al., 1983). 

Meridional smoothing is then applied through a combination of four passes of a 
low pass (2,3,2) filter and one pass of an inverse (-1,5, -1) filter. The 
frequency response of the meridional smoothing for selected wavelengths is 

given in Table 1. 
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Table 1 


Frequency Response 

Wavelength (km) 

12,000 

10,000 

8,000 

6,000 

4.000 

3.000 

2.000 

‘Meridional spatial filter composed o 
filter and one pass of an inverse (- 


f Meridional Filtering 1 

X of original amplitude retained 

96 

95 

92 

85 

69 

49 

16 

four passes of a low pass (2,3,2) 
,5,-1) filter . 


The same meridional smoothing has been applied to the zonally averaged 
distributions with the following exception. Mass -weighted zonal averaging 
acts to smooth data in proportion to the meridional distribution of mass which 
tends to vary proportionally to the cosine of latitude. Recognizing this 
inherent smoothing, the application of the meridional smoothing was made 
proportional to the absolute value of the sine of latitude thereby 
complimenting the smoothing inherent in the mass-weighted average. The 
weighting by the sine of latitude results in an increase of smoothing with 
l a ^itude, from zero at the equator to a maximum at the poles. 

4. Uncertainties in Estimation of Heating 

The accuracy of the heating estimates for the year of the GWE presented 
in this atlas depends on many factors. The ECMWF Level III data are prepared 
through use of an assimilation model. Therefore, the accuracy of the heating 
distributions , being determined from a combination of observations and model 
generated information, depends on the quality of the observations, the 
assimilation model and the diagnostic method utilized. Since a detailed 
analysis of the errors was not feasible, the following discussion is of a 
general nature that aims simply to direct attention to potential sources of 
errors in the analyses. 


In estimating the four -dimensional distribution of heating within 
isentropic coordinates, the interpolation of the Level III data to form an 
isentropic data set and the creation of the surface data introduces 
uncertainty. The lack of information concerning the diabatic mass flux at the 
upper level of the analyses and the uncertainty regarding the adjustment 
employed to remove the bias and insure mass balance are also potential sources 
of errors in the computations. 

In estimating the isobaric heating distributions in this atlas, the 
interpolation of heating from isentropic to isobaric coordinates using time- 
averaged data may not truly reflect the time-averaged isobaric profile of 
heating in some regions due to the transient variation of pressure on 
isentropic surfaces. In transient baroclinic disturbances, latent heating 
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likely occurs at higher pressure (i.e., warmer temperatures) than indicated by 
the time-averaged pressure (i.e., a "mean" temperature) for the isentropic 
level to which it is assigned. 


If the primary objective was to determine the structure of heating 
relative to isobar ic coordinates, optimum estimation suggests that the 
interpolation of the distribution of heating from isentropic to isobaric 
coordinates be carried out for each observational time (as in the more limited 
study in time by Schaack et al . (1990)) rather than interpolating the time- 
averaged isentropic distribution of heating to isobaric coordinates employing 
the time-averaged pressure distribution (as in this atlas). Since the primary 
focus of our research has been to determine the heating distribution relative 
to isentropic coordinates no attempt was made to interpolate from isentropic 
to isobaric coordinates for each time period of the GWE year; a step which 
would have required an extensive reprocessing of the entire GWE data set. 


The differences between results employing interpolation of heating from 
isentropic to isobaric coordinates for each time period (Schaack et al., 1 ) 

versus interpolation of time-averaged structure lie within the realm of 
uncertainty. Schaack et al . show results for corresponding isobaric layers 
for January, April, July and October 1979. The most notable differences 
between estimates occur in January over the oceanic storm tracks of the 
Northern Hemisphere. These regions are characterized by the combination o 
sharp temperature gradients associated with land-ocean contrasts and the 
passage of baroclinic disturbances. In these regions, the results of Schaack 
et al. (1990) generally show stronger heating in the lower and middle 
troposphere while the present results are generally stronger in the upper 
troposphere . 

Johnson et al . (1987) presented layer -averaged results for January and 

July 1979 and zonally averaged results for January, April, July and October 
1979 computed from the ECMWF GWE data set. In the free troposphere the 
methods used in Johnson et al. and the present study are the same. However, 
the present study uses a different method to estimate pJ e (an inverse static 
stability measure as shown in (2)) near the earth's surface. The effect of 
this difference is confined to the low troposphere primarily in high y 
baroclinic regions. The estimates of P J e in these regions in the present 
study tend to be larger than those of Johnson et al . (1987), thereby reducing 
the magnitude of heating (see Eqs . 3 and 4). Comparison of corresponding 
distributions shows the only notable differences between estimates in global 
isobaric layers occur in the surface- to- 800 mb layer of January over the 
western North Atlantic where the results of Johnson et al . (1987) show 
stronger heating. The largest differences between the two studies are in the 
zonally averaged heating distributions where the results of Johnson et al . 
have stronger low-level heating in middle latitudes, particularly in the 
Northern Hemisphere in January. 


While a number of uncertainties are involved, the apparent spatial and 
temporal continuity of the analyses and their consistency with known. features 
of the large-scale circulation suggests that a realistic four-dimensional 
structure of the global distribution of heating during the GWE emerges. 
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Appendix A: Estimation of Heating 

The time -averaged isentropic mass continuity equation is (Johnson, 

d d 

7~(P J 0) + v 0*(pJ0U) + — (pj 00) = 0 . 

dt e ~ 30 

With the hydrostatic relation 
P J 0 “ -g" 1 dp/30 , 


1980) 

(Al) 


(A2) 


the mass distribution is determined by the pressure distribution. 

The time -averaged diabatic mass flux through an isentropic surface 
obtained by indefinite vertical integration of (Al) is 


• 0-p d 

pj00 " ^0 [ 3^ (pJe) + v e*(pJ 0 U)] d0 , (A3) 

where the diabatic mass flux has been assumed to vanish at the top isentropic 
level, 0 T . The boundary value for diabatic mass flux at the earth's surface, 
obtained from a definite integration of (Al) from the earth's surface at 0. to 
the upper surface is given by 


• . ©X ^ 

pj e e| - r [— (pj 0 ) +v.(pj e u)] d© . (A4) 

u s 0 s dt e 


The mass -weighted time -averaged heating is obtained from the diabatic mass 
flux according to 

A t 

• • 

e - pj 0 0 / pj 0 . 


The estimates of heating at lower levels calculated in this manner from 
the assimilated data are unrealistic at some locations. This results 
primarily from biases in the estimation of horizontal mass divergence from the 
assimilated data. In analyses of fields, whether by subjective or objective 
means, bias errors are readily introduced in estimation of divergence (Schmidt 
and Johnson, 1972). These errors stem primarily from discrete irregularly 
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spaced observations which are sparse in relation to the horizontal variation 
of the velocity and the mass pJ e • While a bias error may only be a minor 
fraction of the true horizontal divergence within a particular layer, such 
errors are prevalent in the assimilated data and accumulate in vertical 
integration. Consequently the integrated bias should be removed. For this 
step, a mass -weighted adjustment is applied to the mass tendency and 
divergence subject to the constraint that the sum of the vertically integrated 
mass tendency and mass divergence must be equal to a particular boun ary va ue 
for diabatic mass flux at the earth's surface. This particular boundary val 
must be defined in order to allow the vertical movement of isentropic surfaces 
relative to the earth's surface in accord with changes of potential 
temperature at the earth's surface. 

An independent estimate of the boundary value of diabatic mass flux is 
given by the temporal variation of surface potential temperature, 6 s (A,<t>,t; 
[equal to 0( A ,<J> , z s , t) ] , from the chain rule 


ae, 


“ P J e t— + u s* v e 0 sl 

OtQ ~ 


(A6) 


where A is longitude, <|> is latitude and the subscript s denotes values at the 
earth's surface. The accuracy of this estimate is dependent on the accuracy 
of estimates of the tendency and advection of 9 s (A,<|»,t) at the earth * 
surface. Biases will stem primarily from inadequate resolution of information 
at the earth's surface. The largest biases will be found during winter 
primarily in the regions of strong surface winds and gradients of e s ,sucha 
over the Kuroshio Current and the Gulf Stream just off the east coasts of Asia 
and North America, respectively. Analysis of estimates from t ^ i S ' ' 
however, establish that the global patterns of the surface distribution of 
heating are much more realistic than those from (A4) . 

With the assumption that the lower boundary value of diabatic mass flux 
by (A6) is an unbiased estimate of the "true" value 


• * 

p j e e| - 

®s 


pJ e 0 s (A,(|>, 


t) 


(A7) 


an estimate of the vertically integrated bias error (6) from the difference 
between (A4) and (A7) is given by 


. * 


6(A,<J>,t) - pJQ0| - P^e® I 


(A8) 


®s e s 


The error per unit mass within the atmospheric column is then defined by 


A Q - 

& — 6 / / T P-Je d® 

0 c 


- -g6 / [p(®x) ' P< e s)l 


(A9) 
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Through subtraction of the bias error, the adjusted time-averaged diabatic 
mass flux through an isentropic surface is expressed by 


• 0T* 

(pj e er - f x [ 
0 


d 

- — (pj 0 ) + v(pj 0 u) ] - pj 0 6 de , 
dt e 


- pj 0 0 + &[p(e T ) - p( 0 ) ] , (A10) 


where the bias error is removed in proportion to the mass within an isentropic 
layer . 
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Fig. 4: Meridional cross sections of annual (A) isentropically and (B) 

isobarically zonally averaged heating (10' 1 K day' 1 ), and (C) meridional 
profile of zonally-vertically averaged heating for December 1978 to November 
1979 (K day ) . Contour interval in (A) and (B) is 0.25 K day’ 1 . 
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J: Seasonal layer -averaged heating (K day' 1 ) for the (A) surface-800 mb, 
(B) 800-600 mb, (C) 600-400 mb and (D) 400-200 mb isobaric layers for December 
1978 to February 1979. Contour interval is 0.5 K day* 1 . 
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Fig. 23: Global distributions of seasonal-averaged (A) surface pressure 

(10 1 mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for June 1979 to August 1979. Contour interval is 
30 mb. 
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Fig. 29: Global distributions of seasonal -averaged (A) surface pressure 

(10 1 mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for September 1979 to November 1979. Contour 
interval is 30 mb. 
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Fig. 31: Monthly vertical -averaged heating (K day l ) for December 1978; 

(A) global, (B) Northern and (C) Southern Hemisphere. Contour interval is 
0.5 K day' 1 . 
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Fig. 35: Global distributions of monthly averaged (A) surface pressure 

(10 1 mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for December 1978. Contour interval is 30 mb. 
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Fig. 41: Global distributions of monthly averaged (A) surface pressure 

(10* mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for January 1979. Contour interval is 30 mb. 
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Fig. 47: Global distributions of monthly averaged (A) surface pressure 

(10 1 mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for February 1979. Contour interval is 30 mb. 
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Fig. 53: Global distributions of monthly averaged (A) surface pressure 

(10^ mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for March 1979. Contour interval is 30 mb. 
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Fig. 63: Monthly layer -averaged hej 
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Fig. 67: Monthly vertical -averaged heating (K day x ) for June 1979; 

(A) global, (B) Northern and (C) Southern Hemisphere. Contour interval is 
0.5 K day' 1 . 

































Fie 71' Global distributions of monthly averaged (A) surface pressure 
(10 i mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for June 1979. Contour interval is 30 mb. 




















Fig. 74: Monthly layer -averaged heating (K day* 1 ) for the (A) surface- 800 mb 

(B) 800-600 mb, (C) 600-400 mb and (D) 400-200 mb isobaric layers for July 
1979. Contour interval is 0.5 K day' 1 . 

















Monthl y layer -averaged heating (K day" 1 ) for the (A) surface -290 K 
(B) 290-310 K, (C) 310-330 K and (D) 330-400 K isentropic layers for July 
1979. Contour interval is 0.5 K day" 1 . 
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Fig. 77: Global distributions of monthly averaged (A) surface pressure 

(10 l mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for July 1979. Contour interval is 30 mb. 
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Fig. 80: Monthly layer -averaged heating (K day* 1 ) for the (A) surface-800 mb 

(B) 800-600 mb, (C) 600-400 mb and (D) 400-200 mb isobaric layers for August 
1979. Contour interval is 0.5 K day' 1 . 
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Fig. 82: Meridional cross sections of monthly (A) isentropically and (B) 

isobarically zonally averaged heating ( 10" 1 K day' 1 ), and (C) meridional 
profile of zonally- vertically averaged heating for August 1979 (K day’ 1 ) 
Contour interval in (A) and (B) is 0.25 K day' 1 . 
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Fig. 87: Monthly layer-averaged heating (K day' 1 ) for the (A) surface-290 K, 
(B) 290-310 K, (C) 310-330 K and (D) 330-400 K isentropic layers for Septembe 
1979. Contour interval is 0.5 K day' 1 . 
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Fie 89’ Global distributions of monthly averaged (A) surface pressure 
no 1 mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for September 1979. Contour interval is m 
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Pig. 95 : Global distributions of monthly averaged (A) surface pressure 

(10 1 mb) and (B) pressure difference ( 10 1 mb) between the surface and the 
400 K isentropic level for October 1979. Contour interval is 30 mb. 
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Fig. 101: Global distributions of monthly averaged (A) surface pressure 

( 10 1 mb) and (B) pressure difference (10 1 mb) between the surface and the 
400 K isentropic level for November 1979. Contour interval is 30 mb. 
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